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We study scaling properties of stochastic aggregation processes in one dimension. Numerical 
simulations for both diffusive and ballistic transport show that the mass distribution is characterized 
by two independent nontrivial exponents corresponding to the survival probability of particles and 
monomers. The overall behavior agrees qualitatively with the mean-field theory. This theory also 
provides a useful approximation for the decay exponents, as well as the limiting mass distribution. 



PACS numbers: 05.40.-a, 05.20.Dd, 82.20.Mj 
I. INTRODUCTION 

In the previous study, we introduced a stochastic 
aggregation process involving both active and passive 
clusters |Q. We generalized Smoluchowski's rate equa- 
tions and obtained exact results for several kernels. In 
this study, we apply stochastic aggregation to reaction- 
diffusion, coarsening, and ballistic agglomeration prob- 
lems. Our goal is to examine the range of validity the 
mean-field results, and to determine whether the overall 
scaling behavior extends to low dimensional systems. 

The rate equations approach is mean-field in nature, 
i.e., it is valid only when spatial correlations are absent. 
Formally, it is applicable only in infinite spatial dimen- 
sion, or in the presence of an effective mixing mechanism. 
This mean-field theory should also be asymptotically ex- 
act when the spatial dimension is sufficiently high. In 
low spatial dimensions, however, significant spatial cor- 
relations eventually develop, and the rate equation ap- 
proach does not apply in the long time limit. We there- 
fore focus on one-dimensional systems where spatial cor- 
relations are most pronounced. 

We performed numerical simulations of stochastic ag- 
gregation processes with both diffusive and ballistic par- 
ticle transport. The simulations show that the scal- 
ing behavior suggested by the mean-field theory is in- 
deed generic, as it extends to one-dimensional systems. 
We find that two nontrivial model-dependent exponents 
characterize the survival probabilities of the particles and 
monomers, respectively. Smoluchowski's theory provides 
reasonable estimates for these exponents. 

Additionally, we studied the limiting mass distribution 
of passive clusters. Surprisingly, over a substantial mass 
range, this distribution depends only weakly on the un- 
derlying transport mechanism. Furthermore, mean-field 
theory provides an excellent approximation for the lim- 
iting mass distribution. 

The rest of this paper is organized as follows. The gen- 
eral scaling behavior is outlined in Sec. II. Predictions of 
the mean-field theory are summarized in Sec. III. Nu- 
merical simulations of stochastic aggregation processes 
with diffusive and ballistic transport mechanisms are de- 
scribed in Sees. IV and V, respectively. A discussion of 
the results is presented in Sec. VI. 



II. SCALING PROPERTIES 

Stochastic aggregation involves two types of clusters: 
active and passive JjJ. Initially, the system consists of 
active monomers only. When two active clusters merge, 
the newly-born aggregate remains active with probabil- 
ity p, or becomes passive (i.e., it never aggregates again) 
with probability q = 1 — p. Eventually, all active clus- 
ters are depleted and the system consists of passive clus- 
ters only. This process can be viewed as an aggregation- 
annihilation process since it interpolates between aggre- 
gation (p = 1) and annihilation (p = 0) ||. 

Quantities of interest include A k (t) and P k it) , the dis- 
tributions of active and passive clusters at time t, as 
well as the final distribution of passive clusters, Pfc(oo). 
As shown in 0, two conservation laws underly this 
system. The first is mass conservation, ^fcLAfc(i) + 
Pk(t)} = const. The second conservation law reflects the 
fact that changes in the overall densities are coupled, 
qA(t) + (1 + q)P(t) = const, where A(t) = J2 A k (t) and 
P{t) = J2Pk(t) are the number densities of active and 
passive clusters, respectively. 

Therefore, it is sufficient to study the time evolution of 
the number density and the mass density of active clus- 
ters, A(t) and M(t) = kAk{t), respectively. The latter 
quantity is the survival probability of an active particle, 
i.e., the probability that it still belongs to an active clus- 
ter at time t. Smoluchowski's theory suggests that both 
quantities decay algebraically in the long time limit 

A{t)~t- V , M(t)~t-^. (1) 

As will be shown below, this as well as other scaling prop- 
erties suggested by this theory hold qualitatively even for 
low dimensional stochastic aggregation processes. While 
the decay exponent v is typically robust in that it de- 
pends only on the major characteristics of the process 
such as the spatial dimension or the transport mecha- 
nism, the exponent ip = ip(p) is non-universal as it de- 
pends on the details of the model, i.e., on the probability 
p. In turn, this implies a non-universal growth law for the 
average mass of an active cluster (k) = M/A ~ t vtyl ~^ . 

For the system to follow a scaling behavior, the aver- 
age mass must be the only relevant scale in the long time 
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limit, and conversely, any scale characterizing the initial 
mass distribution must be "erased" eventually. In other 
words, the mass distribution is characterized by a single 
rescaled variable 



A k (t) ~t v ^-^F(kt v 



(V-i) 



(2) 



with the time dependent prefactor fixed by the decay laws 
©■ 

This scaling behavior is similar to that found for deter- 
ministic aggregation-annihilation processes [§-§ and for 
aggregation-annihilation of domains in coarsening pro- 
cesses |j|-|8). These studies suggest that another inde- 
pendent exponent describes the decay of small clusters. 
Specifically, the monomer density decays according to 
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with a model-dependent exponent 5 = 8{p). The 
monomer density decay reflects the small argument be- 
havior of the scaling function F(£) ~ £ CT with 8 — 1 = 
(1 — tp)(l + er). One of our main results is that the mass 
distribution of active clusters is described by a set of non- 
trivial exponents (ip,8). These exponents can be viewed 
as persistence exponents [^[[lO) as they characterize the 
survival probability of an active particle, and an active 
monomer | flT| . 

Several properties of the scaling exponents are gen- 
eral. For instance, the inequalities ip < 1 < 8 hold 
since A\ < ^ Ak < kA-k ■ The two exponents are 
equal ip = 5 = 1 in the annihilation case (p — 0), since 
Ak(t) = A(t)Sk,i- In the aggregation limit (p = 1) the 
mass density of active clusters is conserved and therefore 
ip = 0. 

We now turn to the mass distribution of passive clus- 
ters. The Smoluchowski theory suggests that the same 
scaling form underlies both mass distributions 
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In contrast with the active cluster distribution, the pas- 
sive cluster distribution approaches a nontrivial final dis- 
tribution Pk(oo). Such a time independent final distribu- 
tion is consistent with the above scaling form only when 
the scaling function diverges, F(£) ~ £ r in the limit 
£ -> 0, with 7 = (2 - V)/(l ~ ip)- As a result, the final 
mass distribution of passive clusters decays algebraically 
in the large mass limit 



Pk(oo) 



with 
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At a given time t, this decay is realized for clusters 
whose mass fc does not exceed the characteristic mass 
k* ~ f^i-VO. Note also that < -0 < 1 implies 
2 < 7 < oo. Generally, the mass conservation restricts 
the large mass decay exponent to 7 > 2. Since the ip 
exponent varies between and 1, we see that the entire 
range of acceptable exponents is realized by tuning the 
probability p. 



III. MEAN-FIELD THEORY 



It is well established that spatial correlations can be 
safely neglected only in spatial dimensions larger than 
some upper critical dimension, d > d c fl^| . For example, 
for irreversible aggregation with mass-independent diffu- 
sion and reaction rates, one has d c = 2; for a general ag- 
gregation process, however, the upper critical dimension 
may be arbitrarily large |13j| . Below the upper critical 
dimension, substantial spatial correlations develop, and 
the most important features including the scaling expo- 
nents and the scaling functions are changed. Generally, 
the lower the spatial dimension, the larger the difference 
with the mean-field predictions. 

Although the Smoluchowski rate equations approach 
does not apply in low spatial dimensions, it can still serve 
as a useful approximation after an appropriate modifica- 
tion. This can be accomplished by replacing the overall 
reaction rate with an effective density-dependent reaction 
rate r = r{A) 



dA k 
dt 

dP k 
dt 



r \ \ A * A J - AkA 
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We are primarily interested in situations where aggre- 
gation is independent of the mass, and therefore we 
use a mass independent rate kernel. The reaction rate 
r{A) is model dependent. In reaction-diffusion processes, 
the reaction rate d ecay s algebraically with the density 
(see, e.g., Refs. JlJJilf). Assuming r(A) ~ A a yields 
&rr ~ ~A a+2 , and consequently, the density decay expo- 
nent is found 
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In general, a reduction in the reaction rate, i.e., a > 0, 
leads to a slowing down in the density decay rate, v < 1. 
Apart from the change in v, all other aspects of this 
approximation are identical to the Smoluchowski theory 
with a constant rate kernel. Indeed, the above rate equa- 
tions reduce to the Smoluchowski's rate equations with 
a redefined time variable, t — > t = J Q dt' r{t'). In partic- 
ular, the scaling exponents ip and 5 are independent of 
a: 
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One can verify the expected limiting behaviors ^(1) = 0, 
and ip(0) = S(0) = 1. Furthermore, the scaling functions 
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are as in the constant kernel solution 0]], and for exam- 
ple, F(£) is purely exponential. The corresponding small 
argument exponents 7 = 2/p and a — follow from tj) 
and 5 using the aforementioned scaling relations. The 
final mass distribution of passive clusters is independent 
of the reaction rate r M 
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q T(l + 2/p)T(k) 
p T(k + 2/p) 
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Below, we compare these mean-field predictions with 
simulation results for one-dimensional stochastic aggre- 
gation where spatial correlations are most pronounced. 
We also examine the role of the aggregates' transport 
mechanism by considering both diffusive and ballistic 
transport. 



IV. DIFFUSIVE TRANSPORT 

In diffusive stochastic aggregation, identical particles 
are placed onto a e?-dimensional lattice. All particles per- 
form independent random walks, i.e., they hop to a ran- 
domly chosen nearest-neighbor site with a constant rate. 
If this site is occupied, the two particles coalesce irre- 
versibly, and with probability p the resulting aggregate 
remains active, while with probability q = 1 — p it be- 
comes passive. Effectively, passive particles are removed 
from the system. 

In the case of single-species reaction diffusion pro- 
cesses, the effective reaction rate can be obtained from 
dimensional analysis. Eq. @ implies [r] = [i] d [T] _1 , 
and since the reaction rate can only be a function of 
the diffusion coefficient [D] = [-L] 2 [T] _1 and the den- 
sity [A] = [L]- d , one finds r oc DA&-<W d . Hence, 
a = (2 — d)/d and Eq. (Q) yields the correct decay ex- 
ponents v = d/2 below the upper critical dimension 
4-2. 

To examine the above scaling picture we performed 
numerical simulations of diffusive stochastic aggregation 
processes in one dimension. Unless noted otherwise, the 
data was obtained from an average over 10 independent 
realizations in a system of size L = 10 7 with periodic 
boundary conditions. Initially, all sites were occupied. 
First, we verified that the number density, the mass den- 
sity, and the monomer density indeed decay algebraically 
in the long time limit, in accord with Eqs. (|l|) and 
The case p = 1/2 is shown in Fig. 1, and the corre- 
sponding decay exponents were found: v = 0.500(1), 
tp = 0.6193(3), and S = 1.460(2). Mean-field theory cor- 
rectly predicts v = 1/2. Furthermore, the predictions 
tp = 2/3 and 5 — 4/3 provide a reasonable approxi- 
mation. One can compare with the case of disordered 
(Sinai) diffusion where a real-space decimation proce- 
dure |ll| was used to determine exact values of these 
exponents ||. Remarkably, the disorder case exponent 
tp = 0.61937 is in excellent agreement with the simula- 
tion value. There is a small discrepancy with the second 



exponent 8 — 1.47041. Additionally, we verified that the 
densities of active and passive clusters follow the scaling 
forms of Eqs. (||) and (|4|), respectively (see Fig. 2). In 
agreement with the mean-field theory, the scaling func- 
tions decay exponentially for large masses. 
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Fig.l The number density, the mass density, and the 
monomer density versus time for p = 1/2. 
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Fig. 2 Scaling of the active and passive mass distributions. 
Shown are the scaling functions = £"( 2_ 'W Ak(t), and 

G(£) = fP-^Pkit), versus the scaling variable £ = H' 41 -^ 
at three different times t = 10 4 , 10 s , and 10 6 . Different scales 
correspond to F(£) and £?(£) in the main figure since the lat- 
ter diverges at the origin. The data represent an average over 
10 3 independent realizations in a system of size L = 10 6 for 
the case p — 1/2. The exponent value tj) — 0.619 was used. 
The tail of the distribution is shown in the inset. 

We also studied how the exponents vary with the prob- 
ability p, as shown in Figs. 3 and 4. The exact exponents 
found for the disordered case by Le Doussal and Monthus 
B provide an excellent approximation (within 0.1%) for 
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ip. In the case of S, a different behavior emerges in the 
aggregation limit, p — > 1, where the exact value is S = 3 
p7[ , and the disagreement with both mean-field theory 
and the disordered case are most pronounced. 




Fig. 3 The exponent ip versus p. Monte Carlo simulation re- 
sults for the pure diffusion case are compared with the mean- 
field theory (^) and the exact value for the disordered case. 
The latter is obtained from 2-fo( — g^-, 2ip, 2) = ||, where 
2i ? o(a,fe, z) is the confluent hypergeometric function. 
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Fig. 4 The exponent S versus p. Monte Carlo simulation re- 
sults for the pure diffusion case are compared with the mean- 
field theory ^ and the exact value for the disordered case 
obtained from 2-Fo(-2§5§, 25, 2) = §. 

The above scaling arguments suggest that the limiting 
mass distribution of passive clusters decays algebraically 
with the exponent 7 = (2 — ijj) / (1 — ijj) . For p = 1/2, one 
therefore expects 7 = 3.627 (compare with 7 = 3.62722 
and 7 = 4, predicted by the disordered case and the 
mean- field theory). This corresponds to a very strong 
suppression of large masses, and therefore, it is much 



more difficult to confirm this behavior numerically. Nev- 
ertheless, our simulations (Fig. 5) are consistent with a 
power law decay with an exponent 7 = 3.6. 

In one dimension, the diffusion-controlled stochastic 
aggregation is equivalent to the Potts model with zero- 
temperature Glauber dynamics Jig] . For the Q-state 
Potts model with spatially uncorrelated initial condi- 
tions, aggregation of domain walls occurs with proba- 
bility p — i an d annihilation occurs with probability 

q = Therefore, the above can be reformulated in 

terms of domain walls rather than aggregates. In the 
coarsening context, the domain wall mass (or number) 
distribution is dual to domain number distribution 



10" 
1Cf 1 
10 2 

10- 3 
V 4 
10- 5 

10- 6 

1(T 



\\ 










\\\ 


t=10 4 


\\\ 


t=10 5 


\ \\V 

\ \\ 


t=10 6 


\ \ \ 


slope=-3.6 


'-. \ \ 

\ \ \ 

\ \ \ 



10" 



10' 

k 



10' 



Fig. 5 The distribution Pk(i) versus k for three different times 
t = f0 4 , 10 5 and 10 6 . The typical mass at these three times 
is proportional to k* = oc 6, fO, and f6, respectively. 

Hence, the distribution is fully developed only over a short 
range of masses. The data represents an average over 10 3 
realizations in a system of size L = fO 6 with p= 1/2. 



V. BALLISTIC TRANSPORT 

The situation when particles move ballistically involves 
several complications. First, while the annihilation limit 
is uniquely defined [^9|-|3), the aggregation limit has 
various realizations. In traffic flows, the velocity of a 
newly-born cluster is the smaller of the two velocities 
p4| , while in application to astrophysics and granular 
gases the velocity follows from momentum conservation 
p5|, p6|. Second, numerical results for the annihilation 
case po|] and analytical results for the traffic case |2~i| ] 
show that the initial conditions are remembered forever, 
in contrast with the diffusive case. Specifically, the small 
velocity characteristics of the initial velocity distribution 
influence the long time asymptotic behavior, including 
the scaling exponents. 

We consider the momentum conserving case, also 
known as "ballistic aggregation" or "sticky gas" |27] [33| . 
The initial velocities are assigned according to the dis- 
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tribution Po (v) . The mass (momentum) of a newly-born 
cluster is equal to the sum of masses (momenta) of the 
two colliding clusters. After an agglomeration event, the 
newborn particle remains active with probability p, or 
becomes passive with probability q = 1 — p. 

To apply the Smoluchowski rate equations approach, 
we again use dimensional analysis to calculate the decay 
exponent v. The collision rate is r ~ wa d_1 , where v is 
the typical velocity and a is the typical radius of an ag- 
gregate. A particle of radius a contains of the order a d 
monomers whose initial momenta are uncorrelated. Mo- 
mentum conservation therefore implies v ~ a~ d / 2 . Using 

a d ~ M/A ~ A^- 1 gives the collision rate r ~ a ^ 2 ^ 2 ~ 
A {d-2)^-i)/2d t From Eq (|) one findg 

" = d+2 + t(d~2y (10) 

with ip given by Eq. (^). Apart from the exponent v, 
features such as the exponential mass distribution and 
the exponents ip and 5 are given by the mean-field the- 
ory outlined above. In two dimensions, the collision rate 
does not depend on ip and hence, the asymptotic be- 
havior A ~ t^ 1 agrees with that found for deterministic 
ballistic agglomeration [^7j. For d =/= 2, stochastic and 
deterministic asymptotics differ: stochasticity enhances 
decay of the number density A for d < 2 and weakens 
it for d > 2. A more detailed mean-field theory can be 
carried. It yields a factorizing joint mass-velocity dis- 
tribution, with an exponential mass distribution, and a 
Gaussian velocity distribution [p8| , pT| . 

In the aggregation case, ip — and therefore the cor- 
rect scaling exponent v — 2d/(d + 2) |^7j is recovered 
from Eq. (|l0|). For the annihilation case, however, initial 
conditions are "remembered" forever and therefore the 
above dimensional arguments no longer hold. The pre- 
dicted exponent in the annihilation case is always mean- 
field v = 1, while one-dimensional numerical simulations 
yield an exponent continuously varying from to 1 de- 
pending on the initial velocity distribution Pq(v), e.g., 
v w 0.8 for uniform initial distributions pc|j2^ ]. 

We have simulated the stochastic aggregation process 
on a one-dimensional ring with 10 6 particles. The initial 
velocity distribution was uniform in [—1,1]. We mea- 
sured the scaling exponent ip via the scaling relation 
M ~ A*, rather than directly versus time, since the ex- 
ponent v{p) is not known analytically. We have found 
that the mean-field prediction, ^ = (2 — 2p)/(2 — p), pro- 
vides a reasonable approximation for the exponent ip, as 
shown in Fig. 6. Furthermore, this approximation should 
improve in higher dimensions. 

We compared the mean-field prediction for the mass 
distribution of passive clusters, Eq. (|]), with the numeri- 
cally obtained distributions in both ballistic and diffusive 
cases. Interestingly, the rate equations provide an excel- 
lent approximation for small and moderate masses (see 
Fig. 7). Given that the discrepancy in ij) is maximal for 
the case p = 1/2, one may expect an even better approx- 
imation for other values of p. Noting the strong decay of 



this distribution, the contribution of very large masses is 
extremely small; for example, Pioo(oo) w 2.4 x 10~ 7 for 
p = 1/2. Hence, the most pronounced part of the distri- 
bution is well approximated by the rate equations theory. 
Surprisingly, the transport mechanism does not play an 
important role as far as the final mass distribution of 
passive clusters is concerned. 




Fig. 6 The scaling exponent if)(p) versus p for ballistic aggre- 
gation compared with the mean-field value of Eq. (fel) . 




Fig. 7 The final distribution of passive clusters for the p — 1/2 
stochastic aggregation with diffusive and ballistic transport. 
Also Shown is the mean-field distribution Pk(oo) — 24/[k(k + 
l)(fc + 2)(fc + 3)]. 



VI. DISCUSSION 

We have investigated diffusion- and ballistic-controlled 
stochastic aggregation in one dimension. We have seen 
that the rate equations approach captures the overall 
scaling behavior and additionally it provides reasonable 
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estimates for the decay exponents. In general, the mass 
distribution is characterized by two nontrivial model- 
dependent decay exponents. 

In the diffusion-controlled case, the exponent ip un ~ 
derlying the survival probability of a particle is in excel- 
lent agreement with the exact results from the disordered 
case. In fact, one cannot dismiss the possibility that the 
disordered and the pure values are identical, based on nu- 
merics alone. However, there is an evident discrepancy 
in the exponent S as the disordered case exponent di- 
verges logarithmically in the aggregation limit. Stochas- 
tic aggregation is equivalent to domain coarsening in the 
zero-temperature Potts-Glauber model. The above expo- 
nents (ip,S) characterize the domain wall number distri- 
bution in analogy with (ipD, <5d) for the domain number 
distribution ||. In the latter case as well, exact values 
calculated for the disordered case provide an excellent 
approximation for the domain exponents. In general, the 
particle survival probability exponent ip is robust, while 
the monomer survival probability exponent 5 is very sen- 
sitive to the details of the process. 

In the ballistic-controlled case, we have shown that 
even in the absence of a consistent mean-field theory, 
some characteristics such as the exponent ip are well ap- 
proximated by the rate equations. Understanding of re- 
action processes with an underlying ballistic transport 
remains largely incomplete. The asymptotic behavior is 
highly sensitive to the initial conditions, and the critical 
dimension is apparently infinite. In fact, exact analyti- 
cal results are available mostly in the aggregation limit 
@,||J3|. 

The most intriguing property of the stochastic aggre- 
gation is the profound lack of universality. Indeed, the 
weak dependence on the transport mechanism is in con- 
trast with the strong dependence on the parameter p. For 
example, our numerical results show that the final distri- 
bution of passive clusters is very close in diffusion- and 
ballistic-controlled situations. Another very impressive 
manifestation of this is the excellent agreement between 
the values of the exponent ip(p) m t ne disordered and 
pure cases. 

This research was supported by the DOE (W-7405-ENG- 
36), NSF (DMR9632059), and ARO (DAAH04-96-1- 
0114). 
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